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Abstract. We report on the structural, electronic and vibrational properties 
of the Gai_a;Ina;N alloy using virtual crystal approximation (VGA) from first- 
principles. A band gap bowing parameter of 3.85 eV is obtained with the TB09 
functional. Phonon density of states shifts to lower frequency as the In content 
is increased. However, VGA ignores disorder effect and is therefore unable to 
describe the broadening of the phonon spectra upon alloying. The role of electron- 
phonon interaction in the temperature dependence of the band gap is also studied 
for GaN, InN and their alloy Gai_a,Ina;N. The calculated zero-point motion 
renormalization and the fitted Varshni parameter over the entire composition 
range are discussed. 
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1. Introduction 

The Gai-xIn^N alloy has attracted great interest due 
to its technological importance in optoelectronics and 
electronics. Direct generation of hydrogen by splitting 
water using solar energy has also been reported for the 
Gai_a:In 2 ;N alloy [1]. The band gap of the Gai-^In^N 
alloy can be tuned to cover nearly the entire solar spec¬ 
trum, from 0.8 eV (InN) to 3.5 eV (GaN), by vary¬ 
ing the composition x. Many basic properties of the 
Gai_a,Ina;N alloy are not well known due to the difculty 
in growing high quality samples. The composition- 
dependence of the band gap, also known as the band 
gap bowing, has been extensively studied in recent 
years. However, consensus on the band gap bowing 
of the Gai_a;Ina;N alloy has not been reached. A bow¬ 
ing parameter b of 1.43 eV is reported from an early 
measurement [2], while more recently a considerably 
larger value of 2.8 eV is experimentally observed [3]. 
Results from first-principles calculations range from 
1.1 eV (HSE06)[4], 1.3 eV (LDA-l/2)[5] and 1.5 eV 
(mBJ)[6] to 2.1 eV (LDA-|-C)[7], depending on the 
exchange-correlation functional used. The enhancing 
effect of Indium clustering on the band gap bowing 
has also been examined [7] [8]. Previous first-principles 
calculations rely on either one reasonably sized spe¬ 
cial quansirandom structure (SQS) supercell or a sta¬ 
tistical ensemble of clusters to account for the com¬ 
positional disorder. The virtual crystal approximation 
(VGA), on the other hand, has the main advantage 
of simplicity in modelling the disordered alloy. VGA 
has demonstrated good accuracy for semiconductor al¬ 
loys and ferroelectric perovskite solid solutions[9][10]. 
It offers an efcient way of studying quantities which are 
computationally more demanding such as vibrational 
properties and electron-phonon interactions. Indeed it 
is found that the inclusion of lattice vibrations, which 
is commonly neglected due to its prohibitive computa¬ 
tional cost, could reduce the order-disorder transition 
temperature by ^30% [11]. 

In the recent past, the temperature dependence of the 
band gaps EgiT) of semiconductors have attracted in¬ 
creasing attentions[12][13][14][15][16][17][18][19]. Ex¬ 
perimentally, this temperature dependence has to be 
determined by photoluminescence (PL) or optical ab¬ 
sorption spectroscopy with caution. For example, it 
is at first reported that in InN the PL peak energy in¬ 
creases monotonically with increasing temperature[20]. 
A subsequent study yields normal redshift of optical 


absorption peak with increasing temperature, and at¬ 
tributes the anomalous blueshift of PL peak with in¬ 
creasing temperature to the fact that the PL peak is 
strongly affected by the localized states [21]. For semi¬ 
conductor alloys, theoretical study on the temperature 
dependence of the band gap is rather at its preliminary 
stage. Attempt has been made to fit the experimental 
Eg{x,T) based on the Varshni equation[22]. To our 
best knowledge, first-principles calculations on Eg{T) 
for GaN, InN and Egix^T) for their alloy Gai-^In^jN 
are still lacking. Of particular importance for the first- 
principles band structure calculations is the zero-point 
motion renormalization (ZPR), while the composition 
dependence of the Varshni paramete is crucial for the 
extrapolation of high-temperature Eg(T) from the ex¬ 
perimental low-temperature values. The present study 
aims at throwing some light on the application of VGA 
in the first-principles calculations of the temperature 
dependence of the band gap. 

2. Computational Method 

First-principles calculations are performed using the 
ABINIT program[23] [24] . In VGA, the potential of 
each atom in the alloy is replaced by a composition- 
weighted average of the potentials of its components. 
We employ the VGA provided by the ABINIT 
program. The pseudopotentials are generated using 
the fhi98PP code[25]. The Ga-3d and In-Ad electrons 
are explicitly included in the valence. The electronic 
wave-functions are expanded in a plane-wave basis with 
a kinetic energy cutoff of 50 Hartree. We use a 6 x 6 x 4 
fc-point mesh for Brillouin-zone sampling. Phonons 
are calculated using density-functional perturbation 
theory (DFPT)[26] on a 8 x 8 x 6 g-point mesh. It 
is well known that the band gap calculated from LDA 
is severely underestimated. In the present study we 
apply the recently proposed TB09 functional[27] (a 
modified version of the Becke-Johnson exchange[28] 
part combined with a LDA PW92 correlation[29] part). 
It allows us to improve the band gap prediction with 
a computational cost only slightly heavier than that 
of LDA. The electron-phonon matrix elements are 
calculated on a 8 x 8 x 6 g-point grid in the adiabatic 
rigid-ion approximation[19]. An imaginary shift [id) of 
0.1 eV is used in the perturbation denominator of the 
sum-over-states to avoid numerical instability. 

The most common empirical relation for the variation 
of the band gap with temperature in semiconductors is 
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the so-called Varshni relation[30]: 

Eg=Eo-aTy{T + l3) ( 1 ) 

where Eq is the band gap at 0 K and a, (3 are fitting 
parameters. Although the Varshni relation bears 
no physical resemblance and incorrectly predicts a 
quadratic temperature dependence for T —>■ 0[31], it 
fits the experimental data remarkably well. In present 
study, we continue to use the Varshni relation for its 
significant popularity in applications. 

3. Results and Discussions 

The calculated lattice constants and band gaps of GaN 
and InN with the TB09 functional are summarized in 
Table I. While LDA severely underestimates the band 
gap for GaN and even incorrectly predicts a metallic 
state for InN, the TB09 functional yields band gaps 
close to the experiments. Our VGA calculations show 
that the dependence of a and c on the In content x 
clearly deviate from the Vegard’s law[32]. In Fig. 1 we 
show the calculated and experimental band gaps. To 
study the effect of local environment relaxations (which 
VGA lacks) on the band gaps, we also perform calcula¬ 
tions using the TB09 functional on fully relaxed SQSs. 
Firstly, the TB09 functional corrects for the LDA band 
gaps independently on x, which further confirms the 
applicability of LDA in calculating the band gap bow¬ 
ing parameter for Gai_a;Ina;N alloy [4] [5]. Secondly, our 
VGA calculations show that the band gap bowing at 
low In content is larger than that at high In content. 
Therefore one single composition-independent bowing 
parameter may not be adequate for accurate descrip¬ 
tion of the nonlinear band gap bowing. If we enforce 
the composition-independent bowing, the correspond¬ 
ing bowing parameter reads 3.85 eV by a least-square 
fitting. Thirdly, the band gap bowing of SQS alloy is 
much smaller than that of VGA alloy. The tendency 
accidentally corresponds to the effect of In clustering 
where the band gap bowing ranges from 2.1 eV for the 
uniform case to 3.9 eV for the clustering case[7]. The 
calculated phonon density of states (DOS) is shown 
in Fig. 2. The peak of phonon DOS shifts to lower 
frequency while the shape of phonon DOS remains un¬ 
changed as In content x is increased. Compared with 
32-atom SQS calculations [11], VGA succeeds in cap¬ 
turing the shift of phonon DOS, but fails in describing 
the broadening of the phonon spectra, since the broad¬ 
ening of phonon DOS upon alloying is associated with 
local environment disorder. 

The temperature dependence of the band gap AEg{T) 
can be decomposed into two parts: the electron-phonon 
interaction (EPI) term and the thermal expansion 
(TE) term. While the EPI contribution to the tem¬ 
perature dependence of the band gap of GaN is calcu- 



Figure 1: Band gaps of Gai_a;Ina;N. Zero band gaps 
represent the incorrect metallic states predicted by 
LDA. SQSs are taken from Ref. [11]. Experimental 
bowing parameters are taken from Ref. [2] {b = 1.43 
eV) and Ref. [3] (6 = 2.8 eV). 
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Figure 2: Phonon density of states of Gai_xIna;N. 


lated from first-principles, the TE contribution[31] is 
calculated from —B{^^)TjQ[‘2aa{T') ac{T')]dT'. 

The bulk modulus S is taken from ab initio 
calculatoins[34], while the pressure coefficient of 


dE 

band gap (-g-^)T is taken from low-temperature PL 


dp 


measurement [35]. The linear thermal expansion co¬ 
efficient a{T) over the entire temperature range is 
described by the experimentally determined Debye 
model [36]. Our calculations agree well with the ex¬ 
perimental data[37], as is shown in Pig. 3. The ZPR 
is as large as -0.15 eV. The main contribution to the 
temperature dependence of the band gap comes from 
the EPI term, which is almost three times of the TE 
term. The fitted Varshni parameters (for only the EPI 
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Table 1: Lattice constants and band gaps of GaN and InN calculated with the TB09 functional. Experimental 


data 

are taken from Ref. [33] and are 

shown in parentheses 


a(A) 

c(A) 

Eg{e\) 

GaN 

3.216 (3.189) 

5.239 (5.185) 

2.93 (3.51) 

InN 

3.521 (3.545) 

5.692 (5.703) 

0.57 (0.78) 


term) a and f reads 0.51 meV/K and 745 K respec¬ 
tively. A simple average of the diverse experimental 
data suggests a=0.91 meV/K and /3=830 K[38], while 
a more recent measurement yields a in the range of 
0.54-0.63 meV/K and f in the range of 700-745K. 



Figure 3: Temperature dependence of the band gap 
(referenced to zero) of GaN. Experimental data are 
taken from Ref. [37]. 

The calculated {x,T) dependence of the band gap is 
shown in Fig. 4. Due to the LDA band gap prob¬ 
lem (see Fig.l), we focus on the Ga-rich (x <0.5) con¬ 
tents. The ZPR and the Varshni parameter a decrease 
in magnitude as x is increased because of the larger 
atomic mass of In, indicating a decrease in the strength 
of EPI. We describe the deviation of the composition 
dependence from linearity by a bowing term —6a;(l—a;), 
as is shown in Fig. 5. The bowing parameter for the 
ZPR is —0.1 eV, which should be taken into account in 
the calculated band gap bowing parameter (in present 
study 3.85 eV). The extrapolated ZPR for InN is - 
0.036 eV, considerably smaller than that for GaN. The 
extrapolated a for InN is 0.096 meV/K, while the ex¬ 
perimentally measured value reads 0.414 meV/K[21]. 
The large difference indicates a significant role of TE 
in the Eg{T) of InN, contrary to the situation for GaN 
where EPI dominates over TE. 



Figure 4: Temperature dependence of the band gaps 
of Gai_a;Ina;N alloy at x <0.5. 



Figure 5: Composition dependence of the ZPR and the 
Varshni a. 

4. Conclusions 

In conclusion, we have studied the structural, 
electronic and vibrational properties of the Gai_xIn 2 ,N 
alloy using VGA from first-principles. We obtain a 
band gap bowing parameter of 3.85 eV with the TB09 
functional. VGA succeeds in capturing the shift of 
phonon DOS, but fails in describing the broadening 
of the phonon spectra. We have also studied the role 
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of EPI in the temperature dependence of the band 
gaps for GaN, InN and their alloy Gai_ 3 ;In^N. For 
GaN EPI plays the dominant role, while for InN TE 
could contribute significantly. The calculated ZPR is 
important for modifying theoretical zero-temperature 
band gap, while the fitted Varshni parameter a is 
crucial for extrapolating high-temperatnre band gap 
experimentally. 
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